Temperature dependence of the coercive field in single-domain particle systems 
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The magnetic properties of CU97C03 and CugoCoio granular alloys were measured over a wide 
temperature range (2 to 300 K). The measurements show an unusual temperature dependence of 
the coercive field. A generalized model is proposed and explains well the experimental behavior 
over a wide temperature range. The coexistence of blocked and unblocked particles for a given 
temperature rises difficulties that are solved here by introducing a temperature dependent blocking 
temperature. An empirical factor gamma (7) arise from the model and is directly related to the 
particle interactions. The proposed generalized model describes well the experimental results and 
can be applied to other single-domain particle system. 
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I. INTRODUCTION 

Since Neel'si pioneering work, the magnetic proper- 
ties of single-domain particles have received considerable 
attention,^ with both technological and academic moti- 
vations. A complete understanding of the magnetic prop- 
erties of nanoscopic systems is not simple, in particular 
because of the complexity of real nanoparticle assemblies, 
involving particle size distributions, magnetic interparti- 
cle interactions, and magnetic anisotropy. Despite these 
difficulties, many theoretical and experimental investiga- 
tions of such systems have provided important contribu- 
tions to the understanding of the magnetism in granular 
systems i2 

An important contribution to the understanding of the 
magnetic behavior of nanoparticles was given by Bean 
and Livingston'^ assuming an assembly of noninteracting 
single-domain particles with uniaxial anisotropy. This 
study was based on the Neel relaxation time 
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where the characteristic time constant tq is normally 
taken in the range 10^^^ — 10^^ s, fcs is the Boltzmann 
constant, K the uniaxial anisotropy constant, and V the 
particle volume. KV represents the energy barrier be- 
tween two easy directions. According to Bean and Liv- 
ingston, at a given observation time (robs) there is a crit- 
ical temperature, called the blocking temperature Tb, 
given by 
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above which the magnetization reversal of an assembly 
of identical (same volume) single-domain particles goes 
from blocked (having hysteresis) to superparamagnetic- 
type behavior. 

The time and temperature behavior of the magnetic 
state of such assembly is understood on the basis of Neel 



relaxation and the Bean-Livingston criterion (Eqs. 
and lO). Within this framework the coercive field is 
expected to decay with the square root of temperature, 
reaching zero in thermal equilibrium.'^ As the effects of 
size distribution and interparticle interactions cannot be 
easily included in this model, this thermal dependence 
of He is widely used. The discrepancies between theory 
and experiment are usually attributed to the presence 
of interparticle interactions and a distribution of particle 
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The effect of dipolar interactions on the coercive field 
and remanence of a monodisperse assembly has been 
recently studied by Kechakros and Trohidou^ using a 
Monte Carlo simulation. These authors show that the 
dipolar interaction slows down the decay of the rema- 
nence and coercive field with temperature. The ferro- 
magnetic characteristic (hysteresis) persists at tempera- 
tures higher than the blocking temperature Tb- However, 
the existence of a size distribution in real systems jeopar- 
dizes the observations of such behavior in experimental 
data. 

Studies considering size distribution have been lim- 
ited to the contribution of superparamagnetic particles, 
whose relative fraction increases with temperature. Al- 
though this contribution has been recognized to be an 
important factor affecting the coercive field, a complete 
description of this effect is still open.^'^'^ As an exam- 
ple, let us consider Cu-Co granular systems that presents 
a narrow hysteresis loop, with a coercive field smaller 
than 1000 Oe even at low temperatures. The isother- 
mal magnetization curve suggests a superparamagnetic 
behavior, while the coercive field has an unusual temper- 
ature dependence. This behavior has been qualitatively 
attributed to interactionsiS or size distribution effectsiii 
Nevertheless, this Hc{T) cannot be quantitatively ex- 
plained by the current models in the whole temperature 
range. 

The combined effects of size distribution and interpar- 
ticle interactions on the magnetic behavior is very com- 
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plex, because both affect the energy barriers of the sys- 
tem. Due to this complexity, the size distributions of such 
systems obtained by different mag netic methodsiSiiiJiJi 
include both effects and may be considered as an 'effec- 
tive size distributions'. 

In this work we develop a generalized method to de- 
scribe the temperature dependence of the coercive field 
of single-domain particle systems. We shall focus on 
the specific samples CugyCoa and CugoCoio granular 
ribbons (iSiiiii^ The novelty in this approach consists in 
the use of a proper mean blocking temperature, which 
is temperature dependent. The temperature behavior of 
the coercive field of the studied systems has been suc- 
cessfully described in a wide temperature range and may 
be applied to other systems. 



II. THE GENERALIZED MODEL 

Considering an assembly of noninteracting single- 
domain particles with uniaxial anisotropy and particle 
volume V, the coercive field in the temperature range 
from to Tb, i.e., when all particles remain blocked, fol- 
lows the well-known relation: 



T 



(3) 



Here Ms is the saturation magnetization and a = 1 if 
the particle easy axes are aligned or a = 0.48 if randomly 
oriented i^ii^ 

Some difficulties arise when Hc{T) is calculated on a 
system with several particle sizes: (i) the distribution 
of sizes yields a corresponding distribution of Tb] (m) 
the magnetization of superparamagnetic particles, whose 
relative fraction increases with temperature, makes the 
average coercive field smaller than those for the set of 
particles that remain blocked. These two points will be 
discussed in the following sections. 



Influence of size distribution 



This artifice may lead to a bad agreement with the ex- 
perimental results because it does not take into account 
point {i) conveniently and neglects point (m) of the above 
discussion. Often, a reasonable agreement with experi- 
mental data using this approach is obtained only for T 
well below (Tb),'^- i.e., when there is a great fraction 
of blocked particles. In this work, the assumption (i) 
and (ii) are used to determine a more realistic average 
blocking temperature. 

The idea consist in taking into account only the 
blocked particles at a given temperature (T), i.e., only 
particles with Tb > T. Thus, Eq. ||2Jl should be rewrit- 
ten as 



HcB = a- 
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where {Tb)'y is the average blocking temperature, defined 
by: 



{Tb)t 



j^TBf{TB)dTB 

J^f{TB)dTB ' 
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where the new limits of integration takes into account 
only the blocked particles. 



B. Influence of superparamagnetic particles 

The distribution of particle sizes causes a temperature 
dependence of the coexistence of both (1) blocked and 
(2) superparamagnetic particles. The influence of su- 
perparamagnetic particles on the coercive field was ex- 
plicitly taken into account by Kneller and Luborskyji 
They considered that the magnetization curves of com- 
ponents (1) and (2) are linear for H < Hcb (see Fig.^. 
Hence Mi = Mr + {Mr / Hcb)H for component (1) and 
M2 — XsH for component (2), where Mr is the rema- 
nence, H the applied magnetic field and xs the super- 
paramagnetic susceptibility. These two components may 
be linearly superposed and the average coercive field be- 
comes 



In granular systems there will be a particle size dis- 
tribution which gives rise to a distribution of blocking 
temperatures TB. One can 

The conventional way used to take into account the 
particle size distribution on the coercive field is to define 
a mean blocking temperature {Tb) corresponding to a 
mean volume {Vm) and then consider Hc{T) in Eq. Q 
with Tb replaced by the average blocking temperature of 
all the particlesiiSfliSi 
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Io"°fiTB)dTB ' 
where / {Tb) is the distribution of blocking temperatures 
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C. Determination of {Hc)rp 

In order to obtain {Hc)^- by Eq. {Tj), three terms must 
be evaluated: Mr{T), xs{T), and Hcb{T) all determined 
from experiments. For the last two terms, a good deter- 
mination of /{Tb) is needed. The isothermal remanent 
magnetization is related to /{Tb) according to>i^ 
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FIG. 1: Contribution of the superparamagnetic particles to 
the coercive field. 

Clearly the derivative of Eq. ((HJ is a direct measure of the 
blocking temperature distribution /(Tb) i.e., dIr/dT cx 

The superparamagnetic susceptibility has two contri- 
butions: isolated and groups of a few Co atoms, xsAiT), 
and particles, xsp{T). The later is the initial suscepti- 
bility in the low field Hmit given by MgY/SksT. For a 
system with nonuniform particle sizes, it can be calcu- 
lated as: 

i\^2 fVc{T) 

where the linear relation between Tb and V (Eq. ||2Jl) was 
used to express XSpiT) in terms of Tg, and the critical 
volume Vc above which the particle is blocked. Thus, the 
superparamagnetic susceptibility can be written as 

Xs{T) = 1^ f^TBf{TB)dTB + ^, (10) 

where the second term is the contribution of groups of a 
few Co atoms defined in terms of the Curie constant C . 
Finally Hcb{T) is determined from the Eqs. ^ and ((HJ. 

Alternatively, /{Tb) can also be determined from zero 
field cooled (ZFC) and field cooled (FC) magnetization 
experiments^^ when one is sure that there are no blocked 
particles at the highest measuring temperature, which is 
not always the case. 

III. EXPERIMENTAL 

We have studied two Cu-Co ribbons of nominal compo- 
sition CU97C03 and CugoCoio in the as cast form. These 
samples were prepared by melt spinning as described in 
reference 

The hysteresis loops M{H) of these samples were mea- 
sured in the temperature range from 2 to 300 K and max- 
imum field of 90 kOe. 



The temperature dependence of the coercive field, and 
remanent magnetization were determined from the hys- 
teresis loops. The saturation magnetization was deter- 
mined by extrapolation of M{l/H) for = 0, at 2 K. 

ZFC curves were measured cooling the system in zero 
magnetic field and measuring during warming with an 
external field applied. The FC curves were measured 
during the cooling procedure with field. 

The Ir{T) curve was measured cooling the systems in 
zero magnetic field from room temperature. At each tem- 
perature the samples were submitted to a field of 7 kOe 
which gives a negligible remanent field in the supercon- 
ducting magnet and is well above the field where the hys- 
teresis loop closes. Following the same procedure used by 
Chantrell et al}^, Ir{T) was determined by waiting 100 s 
after the field is set to zero. 

All measurements were performed using a Quantum 
Design Physical Property Measurement System (PPMS) 
model 6000. 



IV. RESULTS AND DISCUSSION 

A. Distribution of energy barriers and system 
nanostructure 

The morphology of a granular system consisting of 
magnetic precipitates plays an important role on the 
macroscopic behavior. Many efforts were made to in- 
vestigate the morphology of granules in Cu-Co alloys. 
It is very difficult to use electron microscopy to identify 
Co granules in a Cu matrix due to the almost identical 
atomic scattering factor of Cu and Co atoms and very 
similar lattice parameter between the Cu matrix and Co 
granules'^^. Other experimental techniques could only be 
used to estimate the average size distribution^^'^^ "^^. 

In many nanoparticle systems the experimental mag- 
netization curves are remarkably well fitted using a super- 
position of properly weighted Langevin curves, usually 
considering a log-normal distributioniSS Making the as- 
sumptions that the magnetization of each spherical par- 
ticle is independent of its volume, the particle size distri- 
bution may be obtained. In fact, the M{H) curves for the 
sample CugoCoio shown in Fig.|2[a) exhibits a superpara- 
magnetic shape, the same occurring also for CU97C03. 

However, M{H) is more sensitive to the average mag- 
netic moment (or average particle size) than the width 
or distribution shapeiS^ For this reason the use M{H) 
to determine /(Tb) is not a good choice. In addition, a 
closer look to the small hysteresis shows that while the 
remanence decays with temperature, the coercive field at 
20 K is smaller than at 4 and 300 K (see Fig.j^Ib)). This 
issue will be discussed in next section. 

We first used ZFC/FC curves (shown in Fig.|3J) as usu- 
ally done to determine the distribution /{Tb). Is is clear 
that the warming and cooling curves do not close as to 
the highest measured temperature. The histeresis in the 
M{H) curve (see fig.|21b) at 300 K show the presence of 



4 



blocked particles. The determination of the distribution 
/(Tb) using this ZFC/FC curves would clearly lead to an 
inaccurate result because it would neglect the remaining 
blocked particles above T — 330 K. We found thus better 
to determine the distribution by the use of eq.|Hland teh 
experimental remanence. 
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FIG. 2: (a) Hysteresis loop for the CugoCoio as cast sample 
at room temperature, (b) A detail of the narrow hysteresis 
at different temperatures. 
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FIG. 3: Zero field cooled and field cooled curves for the 
CugoCoio as cast measured at Hdc=WO Oe up to T = 330 K. 
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The lines in Fig. ^ were obtained fitting the experi- 
mental data to the integrated Eq. (|SJ) using the above 
/(Tb) distribution. The free parameters used in the fit 
were Ui, (Tbi), (T2, {Tb2) and the weighting factor A, all 
shown in Table 1. Since KV oc Tb, /{Tb) represents the 
distribution of energy barriers. 



0.01 

0.1 
0.01 
t1E-3 
1E-4 
1E-5 
1E-6 















o Cu.„Co.„ 






Calculated 






a) 




□ ; 
□ 


^ 1 ^ 

\ — fB-~Q 




1 1— H-l 








b) 







1 



T(K) 



FIG. 4: (a) Remanent magnetization for CU97C03 (square 
symbols) and CugoCoio (circles), (b) Distribution calculated 
using Eqs. and lllll . 

We can see in Fig. 01a) that the agreement between 
the theoretical and experimental curves is good for both 
samples. Figure El^b) shows the distribution of energy 
barriers obtained from the derivative of the isothermal 
remanence decay (Eq. (jHl). The obtained /(Ts) confirm 
the inhomogeneous magnetic nanostructure observed by 
other authors.^4,^5,26,27 



TABLE I: Blocking temperature distribution parameters 
Sample ai < Tsi > (K) o-^ < Tb2 > (K) 



C03CU97 f.2 
CoioCugo 0.7 



8.4 
2.7 



0.6 
f.4 



53 
128 



61 
78 



Figure 0f a) shows the temperature dependence of the 
isothermal remanence of CU97C03 and CugoCoio samples. 
CugoCoio presents two inflexion points suggesting two 
mean blocking temperatures {{Tbi) and (Tb2))- This 
leads us to assume /{Tb) as being the sum of two log- 
normal distributions: 



f{TB) 
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exp 
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Some characteristics of the nanostructure can be in- 
ferred from above results. For CugoCoio there is a large 
number of small particles (small energy barriers) respon- 
sible for the low temperature maximum, and a few large 
particles responsible for second peak in the energy bar- 
rier (see Fig. Elb)). For the less concentrated sample 
(CU97C03), the energy barrier of the two groups of par- 
ticles (small and large) is expected to be smaller than 
the observed in more concentrated sample, owing to the 
corresponding reduction in the particle sizes. However, 
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we observe that (Tbi) is bigger than for the more con- 
centrated sample (see Tabled). Such behavior can be 
related to the higher surface anisotropy as the particle 
size decreases 1^ 



B. 



Coercive field 



The temperature dependence of the coercive field 
Hc{T) is shown in Fig.El While the more diluted sample 
presents the usual decrease with T, the CugoCoio sample 
exhibits an unusual Hc{T) with a sharp decrease up to 
20 K, followed by an increase and a maximum around 
180 K. 

The solid lines were calculated by use of Eq. |(7J with 
/(Tb) obtained previously and adjusting the anisotropy 
constant K in Eqs. lO and pOfl . Some additional con- 
siderations has to be made to obtain a good agreement 
with the data. The straight calculation of Hc{T) gives 
the dashed line shown in Fig. IHl It is clear that a better 
agreement could be obtained by a horizontal shift in this 
plot. It is well known that barrier distributions /(Is) 
are field dependent liSiSS The distribution of energy bar- 
riers obtained from magnetization measurements shows 
a temperature shift, which is a function of cither the ap- 
plied field or the magnetization state of the sample>i2iSS 
Allia et ali^ have proposed that the effect of interpar- 
ticle interactions can be pictured by the use of an addi- 
tional temperature term in the Langevin function. In our 
case the results suggest to use /(tTb) instead of /(Ts), 
were 7 is an empirical parameter that takes into account 
the effect of random interactions. This leads to a better 
agreement, as shown by the dotted line curve. So far we 
did not take into account the Curie term in Eq. IjlOl) . By 
doing so, the agreement with experimental data becomes 
excellent as shown by the solid lines in Figs. [S] and |Sl(for 
both samples). The best obtained parameters are shown 
in Table II. 

TABLE II: Best adjusted parameters 



Sample 


K (erg/cm^) 


C (emu*K/Oe*cm^) 


7 


C03CU97 as cast 


5.0 * 10'' 


2.3 


1.4 


CoioCugo as cast 


3.2 * lO** 


1.4 


2.4 



We can see in the Table II that K is smaller for the 
more concentred sample, which is consistent with the 
relative reduction of the surface anisotropy contribution 
with the size of the nanoparticles.^'^"^^ Note that the 
CosCugv sample presents a higher concentration of iso- 
lated groups of a few Co atoms and smaller interaction 
parameter 7. The expected stronger interaction (and 7) 
in the CoioCugo sample may also imply in a reduction of 
K, as suggested by the random anisotropy modeliiS The 
application of this model to a system with a negligible 
interaction provides a good description of IIc(T) with a 
gamma temperature shift in the distribution equal to 1 
(7 = 1). 
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FIG. 5: Coercive field He vs. temperature for the two sam- 
ples investigated: experimental (symbols) and calculated with 
the generalized model (line). 



The interesting behavior of Hc{T) of the CugoCoio 
sample can be understood in terms of xs{T). Initially 
Hc{T) decreases with temperature due to the unblock- 
ing of the small particles (see Fig. [7|). Then, thermal 
fluctuations leads to a decrease of xs and a consequent 
increase of He with T until the large particles start to 
unblock and He {T) decreases again. In the other sample, 
unblocking occurs more smoothly in whole temperature 
range, due to the relative proximity of the (Tsi) and 
{Tb2), and Hc{T) present the expected decrease with T 
(see Fig. EJa)). 

The above description is particulary adequate for sys- 
tems that present considerable interactions and devia- 
tions of Hc{T) around near (Tb). This was more evident 
in the CugoCoio sample due to its unusual behavior. In 
this respect we present a comparison of different possible 



o 

u 

I 



□ 

□ 


□ Hj, experimetal 
y - 1 and X3 - X sp 


□ ■■ 


y =2.4andx3 = x3p 


: \ □ ■■, 


y = 2.4 and = x^+x^^ 








10 100 




T(K) 



FIG. 6: He calculations and experimental data for CugoCoio. 
The full line is the {Hc}^ obtained considering interaction 
and groups of a few Co atoms. 
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FIG. 7: Calculated xs (dotted line), {Hc)^ (solid line), and 
experimental data (open symbols) vs. temperature of the 
CugoCoio sample. 
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FIG. 8: Coercive field vs. temperature of the CugoCoio sam- 
ple: experimental (open squares), calculated {Hc}rp (solid 
and dotted-dashed lines), and Hcb (dashed and dotted lines). 
For detail, see text. 



and Q , that works for well isolated and narrow distribu- 
tions of sizeS]^ The dashed line is determined by consid- 
ering the temperature dependent average blocking tem- 
perature (Eqs. l|2Jl and and shows clearly the need 
to include the superparamagnetic correction. By con- 
sidering further this correction with Eq. and Q the 
dashed-dotted line is obtained with an excellent agree- 
ment at low temperatures, where most of the particles 
are still blocked. When we take into account Eqs. © 
and O, i.e., considering the temperature dependence of 
the average blocking temperature we obtain finally the 
solid line describing well the results in the whole temper- 
ature range. 



V. 



CONCLUSION 



We present a generalized model for the description of 
the thermal dependence of He of granular magnetic sys- 
tems. With this model we describe successfully the tem- 
perature dependence of the coercive field of Co-Cu sam- 
ples. The contribution of superparamagnetic particles 
and the use of a temperature dependent average block- 
ing temperature was shown to be important to describe 
the coercive field in a wide temperature range. The inter- 
actions among the particles were considered through an 
empirical multiplying factor 7 as well as an effective size 
distribution, we believe that with this procedure most of 
the fine magnetic particle systems may be well described. 



scenarios to explain the experimental data (see Fig. |HJ). 
Without taking into account the contribution (i) of su- 
perparamagnetic susceptibility of unblocked particles we 
obtain the curves indicated by the dotted and dashed 
lines. 

The dotted is the standard and widely used Eqs. (PJ 
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